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ABSTRACT 

A numerical model of a steady state, thin accretion disk with a constant 
effective speed of sound is presented. We demonstrate that 'zero torque' inner 
boundary condition is a reasonable approximation provided that the disk 
thickness, including the thickness of the torquing magnetic fields, is small 
everywhere. It is likely that this conclusion is correct also for non-steady disks, 
as long as the total thickness at the sonic point, H c , is much smaller than 
the radius there, r c ps r ms . The very existence of thin disks is not proved 
or disproved in this work, but such disks are believed to exist for moderate 
accretion rates. Within our model there is a small torque at r ms , which may 
increase disk luminosity by several percent. An important result of our analysis 
is that the physically acceptable steady state solutions in our toy model exist 
only for a < 0.14 (100w s /c) 1/3 . 

A significant torque may be applied to a thin disk if there is a large scale 
magnetic field, like in a modified Blandford-Znajek mechanism. 



Subject headings: black hole physics — accretion disks — magnetic fields 



1. Introduction 

Theory of accretion disks is several decades old. With time ever more sophisticated 
and more diverse models of accretion onto black holes have been introduced. However, 
when it comes to modeling disk spectra, conventional steady state, geometrically thin disk 
models are still used, adopting the classical 'no torque' inner boundary condition at the 
marginally stable orbit at r ms (e.g. Blaes, Hubeny, Agol, & Krolik 2001). Recently, the 
'no torque' condition for geometrically thin disks has been challenged by several authors 
(Krolik 1999, Gammie 1999, Agol & Krolik 2000, to be referred as AGK.) One of us (BP) 
did not agree with their claim, and presented simple arguments why the 'no torque' inner 
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boundary condition is natural if an accretion disk is geometrically thin (Paczyhski 2000), 
but the referee could not be convinced. Thanks to the electronic preprint server BP's paper 
is readily accessible to all interested readers, who can judge its validity. 

In this paper we present a more detailed quantitative analysis of the inner boundary 
condition, refining arguments given by BP, and making them more precise. We find that 
AGK were qualitatively correct: a classical thin steady state disk must have some torque 
at the r ms , but the effect is not as strong as claimed by them. Following a short historical 
review of accretion disk theory in Section 2, we discuss the concept of thin disks, and we 
present equations that describe a disk with a fixed value of the effective speed of sound 
in Sec. 3. In Sec. 4 the solutions of these equations and their topologies are studied 
numerically. In Sec. 5 we describe an analytical method to identify solutions of different 
nature. In Sec. 6 we discuss some of the physical results of our work as well as those of 
Gammie (1999). In particular, we explain why the model and the reasoning proposed by 
Gammie, while bringing up an important issue, were quantitatively incorrect. Finally, Sec. 
7 concludes this paper. 

2. Historical outline 

A Newtonian theory of geometrically thin accretion disks was developed in classical 
papers by Pringle & Rees (1972), Shakura & Sunyaev (1973) and Lynden-Bell & Pringle 
(1974), generalized for a relativistic case by Novikov & Thorne (1973), and reviewed by 
Pringle (1981). A relativistic effect: the bending of light trajectories leading to direct 
illumination of the disk by itself, was first elaborated and calculated by Cunningham (1975, 
1976). This is marginally important for a disk accreting onto a Schwarzschild black hole, 
but it is important in the Kerr case, introducing significant non-local effects to the energy 
balance at small radii. In these early papers a geometrically thin disk accreting onto a black 
hole had the inner boundary at the marginally stable orbit, r ms , and the gas was freely 
falling on a tight spiral inwards of r ms . 

There was some confusion about the nature of gas flow near the inner disk edge 
(Stoeger 1976). The issue was first understood as a transition from sub-sonic to supersonic 
flow for geometrically thick disks (Abramowicz et al. 1978), and later for geometrically thin 
disks (Muchotrzeb & Paczyhski 1982). Small a-parameter used in these papers resulted in 
saddle-type critical points. Matsumoto et al. (1984) found that for moderately large values 
of a the critical point was nodal-type. These papers, and many more, were reviewed by 
Abramowicz & Kato (1989), who presented a general discussion of the transonic flows in 
accretion disks. In particular, they pointed out that while such a transition is necessary 
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for the existence of a steady state accretion, it does not guarantee that a sensible global 
solution exists (cf. their Fig. 2). 

All the papers written in 1970s and 1980s adopted some form of 'alpha' viscosity, as 
there was no quantitative physical model for the effective transport of angular momentum 
and the dissipation of energy within accretion disks. Perhaps the best qualitative description 
of these processes was given by Galeyev et al. (1979). A major breakthrough was made 
by Balbus & Hawley (1991) and Hawley & Balbus (1991), who rediscovered a powerful 
magneto-rotational instability in weakly magnetized disks, and pointed out its relevance to 
accretion. A flood of papers followed, with powerful computers making it possible to model 
time dependent accretion flows in 2-D (Armitage 1998, Stone et al. 1999, Agol et al. 2001, 
Stone & Pringle 2001, and references therein) and in 3-D (Flemming et al. 2000, Hawley et 
al. 2001, Armitage 2001, Armitage et al. 2001, Reynolds & Armitage 2001, Reynolds et al. 
2001, Sano & Inutsuka 2001, Krolik & Hawley 2002, Hawley & Krolik 2002, Igumenshchev 
et al. 2003, and references therein). 

A lot of semi-analytical and numerical work was done for the ADAF and CDAF models 
(cf. Ball, Narayan & Quataert 2001, and references therein). It appears that at very high 
and at very low accretion rates disks are geometrically thick, while thin disks may exist only 
for moderate accretion rates. If this view is correct then the early work on geometrically 
thin accretion disks may be relevant for 0.01 < L/L Edd < 0.1. 

A very interesting model of a geometrically thin disk with a geometrically thick 
magnetic corona was recently presented by Merloni (2003). Recent reviews of various modes 
of disk accretion onto black holes were provided by Merloni (2002) and by Blaes (2002). 

Yet another way to modify the picture was proposed by Blandford and Znajek (1977), 
who pointed out that a large scale magnetic field may thread a black hole, and it may 
extract its rotational energy. In recent years the original picture was modified, making it 
more likely that the energy extracted from a spinning black hole is transferred to the inner 
parts of the disk, rather than directly to a distant load (Agol & Krolik 2000, Li 2000a,b, 
Wang et al. 2002, and references therein). Recent XMM observations of an AGN named 
MCG-6-30-15 were claimed to support this possibility (Wilms et al. 2001, Merloni & Fabian 
2003). Cao & Xu (2002) investigated local structure of the accretion flow near the sonic 
point assuming there is a steady torque exerted at the inner disk edge. 
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3. A Thin Disk 

The modern 2-D and 3-D numerical simulations provided the first quantitative and 
meaningful insight into the actual physics of angular momentum transport in accretion 
disks. As computer power increases and the codes become more sophisticated a steady 
progress is to be expected. However, at this time there are significant limitations. The 
numerical accretion flows are relatively thick, the outer boundary is not very far out, the 
cooling processes are not included, the time integrations can be carried out for a rather 
modest multiple of the dynamical time scale. Hence, when it comes to modeling disk 
spectra, conventional steady state, geometrically thin disk models are still commonly used 
adopting the classical 'no torque' inner boundary condition at r ms (e.g. Blaes, Hubeny, Agol, 
& Krolik 2001). The current ADAF or CDAF or variants of Blandford - Znajek models 
do not exceed a toy model level. This is both, natural and useful, and the quantitative 
understanding will be improved with time. For the same reason it is useful to have a good 
quantitative understanding of semi-analytical models of geometrically thin accretion disks. 
The simple toy models will remain useful even when the future 3-D numerical approach 
provides full quantitative comparison with the future observations. 

At the toy model level one makes a simplifying assumption in order to construct 
geometrically thin disk, to bypass the issue of thermal (or in general internal) energy 
balance. Krolik (1999) and Gammie (1999) assumed that the relative disk thickness is 
constant and small, H/r = const <C 1, even inwards of the marginally stable orbit, in the 
so called 'plunging region'. We use an alternative assumption: the effective speed of sound 
is constant and small: v s /c = const <C 1. The two assumptions are roughly equivalent as 
a small effective speed of sound implies small disk thickness, and vice versa, as long as the 
flow is in a hydrostatic equilibrium in the 'vertical', i.e. 'z' direction. This holds even in the 
supersonic flow, as long as radial velocity is much smaller than rotation velocity, i.e. as long 
as the accretion time scale is much longer than dynamical time scale. 

In our view the essence of a thin disk concept is the assumption that it is OK to 
integrate disk structure over its thickness, and to consider all important variables to be 
functions of radius only. In particular, thin disk radiates away energy dissipated locally, 
and no advection is allowed. In effect a diversity of local physical processes can be ignored, 
and the conservation laws of mass, angular momentum and energy determine everything. 
A simplifying ad hoc assumption, be it H/r = const < 1, or v c /c = const <C 1, is added in 
order to follow a transition from a subsonic flow at large radii to supersonic infall at small 
radii. 

The argument presented by Paczyhski (2000) in favor of applicability of a no torque 
inner boundary condition was very simple. Geometrically thin, steady state disk accretion 
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onto a black hole was sub-sonic at r > r in , and supersonic at r < r in , with the inner disk 
radius r in located near the marginally stable orbit, at r rns . The specific angular momentum 
at r in is l in , and it satisfies the equation 

— -y— — ~ ct— — — 1, if a«l and ^«1 (1) 

'in Tin Tin 

(cf. eq. 4 in Paczyhski 2000), where l m l in Z ms is the angular momentum integration 
constant, H is the disk thickness, and a is the viscosity parameter. Eq. (1) follows from 
angular momentum conservation in a steady state accretion flow. Current 3-D numerical 
calculations for thick disks and tori indicate a ~ 0.1, with large fluctuations of all physical 
quantities. The inequality — lo)/kn <C 1 is equivalent to the so called 'no torque' inner 
boundary condition, or more precisely to a very small torque at the inner disk edge. This 
conclusion remains valid even for a ~ 1, as long as H in /r in <C 1, i.e. as long as the disk 
remains geometrically thin at the critical point. Note that eq. (1) is local, i.e. it does not 
matter if the parameter a is constant throughout the flow, or does it vary with radius. The 
two essential assumptions are that the disk is geometrically thin, and that the accretion is 
steady state. 

Please note: we do not use eq. (1) in our numerical model calculation, in which we 
make no assumption about a relation between l in , l rns and Iq. The value of the angular 
momentum constant is calculated adopting two assumptions: the disk is 'Keplerian' at large 
radii, and the flow passes smoothly through the effective sonic point. These two conditions 
select a unique value of l , which turns out to be very close to l ms . Also, the effective sonic 
radius r c is found to be close to the marginally stable orbit, r ms . 

What can be wrong with eq. (1)? There are several possibilities. Thin disks, with 
Hin/r in <C 1 may be physically impossible. We do not know if this is true or not, as we 
have neither theoretical nor observational proof either way. Steady state accretion may 
not exist in nature. Current 3-D simulations have indicated that strict steady state flow 
is not possible, but there is a possibility that when averaged over moderate time scale a 
quasi steady state may still be a sensible approximation. The eq. (1) may hold as long as 
the fluctuating disk thickness remains small at all time. This domain is not accessible to 
numerical calculations so far. Next, there is a global problem pointed out by Abramowicz 
& Kato (1989): an accretion flow may pass through a sonic point near r ms but it may not 
continue all the way into a black hole, as shown schematically in their Fig. 2. In fact, we 
could not find any reference from 1980s demonstrating that the supersonic infall for r < r ms 
continues all the way. This was assumed to be obvious and the global flow pattern was not 
verified numerically. 



There is an issue of terminology: what does it mean that a disk is geometrically thin? 
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In the eq. (1) the disk thickness H refers to the structure which carries stresses. If the 
stresses which make the accretion possible are magnetic, as seems to be very likely, then 
H must refer to the geometrical thickness of the magnetic field structure, not just the 
gas layer. The notion that the magnetic thickness may be larger than gas thickness dates 
back at least to Galeyev et al. (1979), and has found some support in recent 2-D and 3-D 
simulations. Even more importantly, all Blandford - Znajek type models are based on large 
scale magnetic structures to transfer momentum and energy. The distinction between the 
magnetic and gas disk thickness was not considered at all in the classical papers about disk 
accretion, and it was not mentioned by Paczyhski (2000). Yet, it is the scale height of the 
stresses that is important for eq. (1). Therefore, we emphasize that the term: 'geometrically 
thin disk' refers to the thickness of magnetic structures responsible for the momentum 
transfer. 

Another source of ambiguity comes with the terms 'pressure' and a. There are many 
different ways in which the parameter a has been defined in the literature. In the original 
work of Shakura & Sunyaev (1973), a was defined as the ratio of tangential stress to the 
pressure of the viscous fluid. Balbus & Hawley (1991), recognized that this stress is in 
fact dominated by the magnetic terms. Early simulations seemed to indicate that the gas 
pressure dominates the magnetic pressure, at least outside the marginally stable orbit. 
Recent 3-D simulations (e.g. Hawley & Krolik 2002) show that magnetic pressure may 
dominate the gas pressure inside the marginally stable orbit. Therefore, we define a in 
terms of total pressure, including magnetic, so it cannot exceed ~ 1. It is not essential for 
our model to have constant value of a, this is just the simplest assumption. 

Another ambiguity is caused by the 'speed of sound'. It is conventionally defined as 
the speed at which gas pressure disturbances propagate, or it may also include the effects 
of radiation pressure. If the magnetic pressure dominates then the Alfven speed becomes 
much larger than the speed of sound. Rather than study the complicated effects of different 
speeds at which different disturbances may propagate we introduce the 'effective speed of 
sound' v s = (P/p) 1 / 2 , where P is the total pressure, which is relevant for the v s . With 
highly tangled magnetic fields this seems to be a reasonable, though perhaps non-orthodox 
definition. It is the effective speed of sound (the fast magnetosonic speed if magnetic 
pressure dominates) that defines the transition from a sub-sonic to a supersonic flow. It is 
beyond the scope of this paper to provide a rigorous justification of our approach, but we 
think this simplification is sensible. For a disk to remain geometrically thin it is essential 
that the effective speed of sound is much smaller than the speed of light, but it does not 
have to be constant. A theory which would allow us to calculate the effective speed of 
sound from the first principles does not exist. Therefore, we consider the simplest possible 
disk model, and we assume that throughout the disk the 'effective speed of sound' remains 
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constant, i.e. v s = const. This assumption is not essential, it is just the simplest. 

We should note that a key ingredient in the above arguments, and in the analysis 
of the following sections, is the assumption of a hydrostatic equilibrium in the direction 
perpendicular to the thin disk. It holds where the accretion time scale is longer than 
dynamical time scale, i.e. where radial velocity is much smaller than rotation velocity. 
The 'vertical' hydrostatic equilibrium holds even while the radial flow is supersonic, 
approximately as long as r > 2r g . The consequence of 'vertical' hydrostatic equilibrium is a 
simple relation between the effective speed of sound v s , rotation velocity v ro t, disk thickness 
H and the radius r: v s /v rot ~ H/r. It follows that if the disk is assumed to be thin, i.e. 
H/r <C 1, then the effective speed of sound must be much smaller than the speed of light, 
i.e. v s /c <C 1. Therefore, the assumption that the disk is thin implies that magnetic energy 
is assumed to be efficiently dissipated. 

Following all these definitions we present a simple thin accretion disk model in order 
to demonstrate that there is no problem of the type envisioned by Abramowicz & Kato 
(1989), i.e. that once the infall becomes supersonic it never 'turns back'. Our simple 
numerical model reproduces the basic features of classical thin disks, and demonstrates that 
the transonic solution may remain unique not only for a saddle, but also for a nodal critical 
point. But first we recall some concepts of the thin disk models of several decades ago. For 
the derivation of all equations the reader may consult the paper by Abramowicz & Kato 
(1989). 

The single most important difference between Newtonian accretion and an accretion 
onto a black hole is the presence of a marginally stable orbit in the latter. Newtonian gravity 
varies strictly as a power of radius, therefore binding energy and the 'Keplerian' angular 
momentum vary as power laws of radius, and this makes self-similarity an acceptable 
simplifying assumption for some accretion flows. This is not the case in general relativity: 
at small radii 'Keplerian' quantities are no longer power laws of radius, they do not even 
vary monotonically. The 'Keplerian' angular momentum and the corresponding binding 
energy reach a minimum at r ms , the radius of a marginally stable orbit. This feature is 
reproduced with a pseudo-Newtonian potential (Paczyhski & Wiita 1980): 

GM 2GM 
r — r g c 2 

For a test particle on a 'Keplerian' circular orbit we have a relation 

cM = GM = v\_ 
dr ( r - r g f r 

where vk is the 'Keplerian' rotational velocity. It follows that 'Keplerian' angular 
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momentum, l K , angular velocity, Qk, and binding energy, e K , are given as 



Ik = v K r = (GMr) 1/2 , Qk — — ( ) ( — ) ■ (4a) 

\r — r g J r \ r A J \r — r g J 

T v\ GM(r-2r g ) 

e K = * + ^ = - —f -f- (46) 

2 2 (r - r g ) 

At very large radii, r/r g ^> 1, the eqs. (2-4) asymptotically become Newtonian, but the 
differences are large at small radii. In particular, a minimum value of 'Keplerian' angular 
momentum and binding energy is reached at r = r ms = 3r g , just as it does in the relativistic 
Schwarzschild case. The orbits with r > r ms are stable, while those with r < r ms are 
unstable, with a marginal stability at r ms . 

Let us consider the simplest disk model with zero pressure. In this case all streamlines 
are identical with particle trajectories, i.e. they are nested 'Keplerian' orbits. This 
structure can be extended to small radii, but all orbits inwards of r ms are unstable: a 
small perturbation makes a particle spiral into a black hole. Therefore, a zero pressure 
disk is truncated at r ms , it cannot extend inwards of r ms . Such a disk has zero geometrical 
thickness, and it is static, i.e. there is no accretion. 

To make disk accretion possible it is necessary to introduce some means of angular 
momentum exchange between nearby streamlines. In general, this will require a non-zero 
total pressure P, a non-zero effective speed of sound v s ~ (P/p) 1 ^ 2 , and a finite disk 
thickness H . For simplicity we assume the effective speed of sound to be constant: 

H v s v 2 P r 

— w — « 1, v s = const, -f- w — - <1, S = / pdz, (5a) 
r vk c 2 pc 2 J 

We remind the reader that pressure P is the total pressure: 

P = Ptot = Pgas + Prad + P-mag + Pturb- (56) 

The effective speed of sound v s becomes close to Alfven speed when the magnetic pressure 
dominates. 

We assume that the stresses responsible for angular momentum transport are 
proportional to the total pressure and the shear rate. The torque exerted by these stresses 
across radius r is 



/GM\ 1 / 2 



The dimensionless factor a is assumed to be less than unity, and for simplicity we take 
it to be constant. Again, just as it was the case with the 'effective speed of sound', the 
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assumption that a is constant is not essential, it is just the simplest. Note, that Vt varies 
monotonically with radius, i.e. the inner disk rotates faster, and angular momentum is 
transported outwards. 

The aP term, no matter how small, redistributes angular momentum within the disk, 
and accretion becomes possible. In a steady state each disk element gradually loses its 
angular momentum, and slowly spirals inwards. For r > r ms the process can be envisioned 
as a motion along ever smaller nearly 'Keplerian' orbits, with gradually decreasing angular 
momentum. This brings matter close to r ms . Following any additional loss of angular 
momentum the gas cannot find any 'Keplerian' orbit and must plunge into the black hole 
along a spiral, approximately conserving angular momentum. 

Let us assume a steady state, with the mass accretion rate M = const < 0, and the 
radial flow velocity v r < assumed to be a function of radius only. The conservation laws 
of mass and angular momentum are 

M = 27rr£t; r , g = (-M) (I - l ) , (7) 

where Iq is the integration constant (cf. Abramowicz & Kato 1989). Combining eqs. (6) 
and (7) we obtain 

(-v r )(l-l ) = arv 2 s -— - — . (8) 



d\nr ) \£Lk- 

It is reasonable to expect that Zo ~ L (to be verified later), and the last equation may be 
used to estimate the radial flow velocity in the nearly 'Keplerian' thin disk for r > r ms and 

\v r \/v s < I. 

Note, that in eq. (8) only two quantities vary a lot within the flow: v r and (/ — lo), 
while the others are either constant or vary slowly. In particular, according to the classical 
theory of thin accretion disks, near the r ms and inwards of r ms angular momentum is almost 
conserved. This implies that 

(-^— )~2, therefore (-v r ) (I - Z ) « 2arv 2 s ( . (9) 
\ amr J \iIk/ 

The numerical solutions of the following sections will verify this assumption. 

We expect that when matter reaches r r ms any additional loss of angular momentum 
puts it on a spiral infall toward the black hole, during which energy and angular momentum 
are approximately conserved, i.e. 

GM . vl . vL c 2 GM 

t = v rot r w l ms = 3V 1.5 , (10) 



r-r g 2 2 16 
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where v rot is the rotational velocity component of the infall. This is verified with the 
numerical integrations in the following sections. 

To calculate a transition from the sub-sonic radial flow to the supersonic flow we need 
the equation of motion (cf. Abramowicz & Kato 1989) 



v 2 dZ dlnQ 



K 



dr ^ ' S dr s dr 

Differentiating the mass conservation law (eq. 7) we obtain 



= 0. 



(11) 



vjd^ = _vl 
E dr r 

and combining eqs. (11) and (12) we find 

d\av. 



dm iv 
dr 



(12) 



2 2 
V~ 



dlnr 



dlnQx 
dlnr 



- 1 



0. 



(13) 



We introduce dimensionless constants and variables defined as 



a = — < 1, 

c 



b = 



lo 



^>o, 

c 



< a < 1, 



r 

x = — . 

r„ 



(14a) 



(146) 



The equation of motion (13) and conservation of angular momentum (9) may be written in 
dimensionless form as 

din v 1 



din a; (v 2 — a 2 ) 



^ - 



x 



2(x-l) 2 
1.5-v/3 b(x - l)v 



+ 2.5 a' 



x - 0.6 
x- 1 



x x l / 2 v — 2\f2 aa 2 (x — 1) 



(14) 
(15) 



4. Critical Points, Physical Solutions and their Topology 

We seek a physical solution which has the properties 

t;/a<l, a; « 1, for i>3, 
t> = a, for x = x c , 



(16a) 
(166) 



v/a > 1, 



for x <C x c . 



(16c) 



i.e. the flow is sub-sonic and the disk is 'Keplerian' at large radii, the flow passes through a 
critical point, and it becomes supersonic at small radii. 

To better understand the behavior of the physical solutions of eq. (14), and to be able 
to integrate it numerically near the critical point, we decompose it into two equations 



which, together with eq. (15), form a two-dimensional autonomous dynamical system. The 
variable t, is a dummy variable and should not be confused with physical time. 

As mentioned earlier, the physical solution must pass a critical point to reach the 
supersonic regime, otherwise the solution would not be single-valued as the right hand side 
of eq. (14) diverges. Now, the critical points of eq. (14) are the fixed points of eqs. (17) 
and an understanding of the nature of these fixed points is necessary for constructing a 
physical solution. 

At any fixed point, the right hand sides of eqs. (17) vanish. Eq. (17b) implies that 
fixed points are always at v — a, as required by eq. (16b). The right hand side of eq. (17a) 
at the critical point, together with eq. (15) (with v — a) lead to a complicated algebraic 
equation for a; at a given b, or vice versa. The analysis is significantly simplified for small 
values of a and the corresponding relation between b and x c is shown in Fig. 1. We see 
that, in general, there are either no fixed points or two fixed points in the system, while at 
the critical value of b = 1 + 0(aa), there is only a single degenerate fixed point at x ~ 3. 

The nature of the fixed points can be investigated analytically by linearizing eqs. (17) 
around each of them. It turns out that the left fixed point (smaller value of x) is always 
saddle-like, i.e. the two eigenvalues have opposite signs. The right fixed point is of the spiral 
type (two complex eigenvalues) when the two fixed points are far apart, but it switches to 
the nodal type (two eigenvalues of the same sign) as the two fixed points get closer. This 
behavior is shown in Figs. 2-5 which is the typical behavior of a saddle-node bifurcation 
(Guckenheimer & Holmes 1983). 

The integrations presented in Figs. 2-5 were calculated numerically, except near the 
critical points, where an analytic expansion was used. All four figures have the same disk 
parameters: a = v s /c = 0.01 and a = 0.15. They differ in the adopted value of the angular 
momentum constant b. The family of solid lines shown in each figure represents possible 




(17a) 



dlnx 



2 2 

— v — a , 




dt 
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solutions of the eq. (14). The solutions of special interest for us are those which satisfy 
boundary conditions (16a,b,c). As it turns out none of these figures has a desired solution, 
but the figures present the changes in the structure of the critical points. 

Fig. 2 presents a saddle critical point at x c = 2.7, and a spiral point at x c = 3.4. The 
Phase Portraits for b > 0.9997 (including b > 1) is qualitatively similar to Fig. 2. The 
Phase Portrait changes when the angular momentum constant b is reduced, as shown in 
Fig. 3, where the two critical points are saddle at x c = 2.94, and nodal at x c = 3.1. When 
b is reduced down to b = 0.9947 the two critical points merge, as shown in Fig. 4. This 
corresponds to the point marked (x c , 1) in Fig. 1. When b is reduced even more then there 
is no critical point, as demonstrated with Fig. 5, and as was anticipated (cf. Fig. 1). 

A prominent feature of Figs. 2-5 is the 'Keplerian' solution. In the sub-sonic regime, 
(w < a), the the right hand side of eq. (14) is very large for small values of a unless uo ~ 1 
(almost 'Keplerian' angular velocity). This implies that the only slowly varying solution 
of eq. (14) is an almost 'Keplerian' one. The radial velocity for a 'Keplerian' solution is 
obtained by setting uj = 1 in eq. (15), 

2^/2aa 2 x 4n . . 

■» = svv(s-i)-iwa +(V >- {) 

Figs. 2-5 show that all the other solutions tend to merge with the 'Keplerian' solution as 
x decreases. This implies that the behavior of the physical solution close to r ms is almost 
independent of the boundary condition at large radii, since all different solutions merge as 
they move inwards. Therefore, it is easy to satisfy the boundary condition given with eq. 
(16a). However, none of the solutions shown in Figs. 2-5 satisfies all three conditions given 
by eqs. (16a,b,c). 

It is clear that a physical solution may pass through a saddle or a nodal fixed point, 
but not a spiral. Furthermore, we assume that a physical solution is analytic. This singles 
out only two ways of passing a nodal fixed point (slow or fast modes), because a linear 
superposition of fast and slow modes is non-analytic (when v is expressed as a function of x). 
In the case of a saddle fixed point, only two solutions (in the direction of two eigenvectors) 
can pass the fixed point. A conclusion that can be drawn from this consideration is that, 
in general, the 'Keplerian' solution does not analytically pass the fixed points and only for 
discrete values of the angular momentum constant b a physical solution is possible. 

Now, we are ready to analyze different topologies that can occur in the parameter 
space of our model. For a given value of a the nature of critical points can be calculated 
analytically for any pair of values of a and x c , or equivalently a and b. Fig. 6 shows the 
possibilities for a = v s /c = 0.01. The horizontal axis, x, is the position of a fixed point and 
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the vertical axis is the parameter a that appears in the viscosity model (eq. 5). The fixed 
points are of the saddle-type, on the left side of the dashed region, of the nodal type inside 
the dashed region, and of the spiral type on the right side of the dashed region. The solid 
curve aBCd shows the position of physical solutions that allow passage of the 'Keplerian' 
curve through the critical point. The critical points present in Figs. 2-4 are indicated in 
Fig. 6 with diamonds numbered 2, 3, 4, 3, and 2, correspondingly, with a — 0.15 for all of 
them. The diamonds with numbers 7, 9, and 11, correspond to critical points presented in 
Figs. 7, 9, 11; they are all located on the solid line aBCd. 

We note that for a smaller than some critical value asN, the physical critical point is 
of the saddle type (segment aB), while for larger values of a it is of the nodal type (segment 
Bd). The value of asN depends on a. For example for a = 0.01 we have cxsn = 0.08 (cf. 
Fig. 6). 

A critical point of the saddle type, shown in Figs. 7 and 8, is the type that is usually 
encountered in astrophysics. The best known examples are the solar wind (Parker 1958) 
and Roche lobe overflow in binary stars. As argued above, when the parameters a and a 
are fixed then physical solution exists only for a unique value of b and the corresponding 
value of x c . 

In the case of a nodal critical point, the passage is possible through fast or slow 1 
directions. However, a generic solution, which is a combination of fast and slow modes, 
passes the critical point in the slow direction since the slow solution dominates the fast one 
close to the critical point. In this sense, a physical solution that passes a nodal point in the 
fast direction, is unique, similar to the saddle critical point. This kind of passage is seen in 
the BC section of the aBCd curve in Fig. 6 and and example of the phase portrait is shown 
in Figs. 9 and 10. Note that, again similar to the saddle point solutions, there is a smooth 
transition from the sub-sonic to the supersonic regime. 

Finally, the physical solution may pass the critical point in the slow direction. However, 
as argued above, this does not fix the value of b since it is the generic behavior of the 
solutions, and only by requiring the analyticity of the physical solution one may determine 
b uniquely. The points on the Cd segment in Fig. 6 are of this type, and an example of the 
phase portraits is plotted in Figs. 11 and 12. 

The transition from the former type of nodal passage to the latter occurs at a^^. For 
example, for a = 0.01, &nn — 0.14 (point C in Fig. 6). The curve aBCd, at the transition 
point C, is tangent to the boundary of the nodal and saddle regions, where the fast and 



The fast/slow direction refers to the direction of the eigenvector with the larger/smaller eigenvalue. 
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slow directions merge. 

The main difference between the last type of passage and the previous two types is 
that there is a sharp change in slope, right after passing the critical point. The reason 
is that the 'Keplerian' curve, which is connected to the slow direction, turns from stable 
(as x decreases) in the sub-sonic regime to unstable in the supersonic regime, and so the 
physical solution departs from the 'Keplerian' curve after the passage. Since the slopes are 
significantly larger far from the 'Keplerian' curve, there is a sharp change in slope. 

While we have described some properties of solutions with a nodal critical point, a 
general discussion of these matters is beyond the scope of this paper. The central issue 
for us is the inner boundary condition for a steady state geometrically thin accretion disk. 
For a given value of the effective speed of sound v s = ac, and for a < 0.14 we found a 
unique numerical solution, in agreement with the finding of Artemova et al. (2001). We 
verified the uniqueness of numerical solution in two ways. First, we began with a guess of 
the location of the critical point x c , and we searched for the value which provided a solution 
close to 'Keplerian' at large radii. This provided a unique value of the angular momentum 
constant b for the given a and a. We also started numerical integrations of the eqs. (17 and 
15) at large radii, where we adopted a 'Keplerian' model. We searched for the value of the 
angular momentum constant b for which the solution would pass through the critical point. 
There was a unique value of b and x c found for the given a and a and it was identical to 
that found with the first method. There was no difference in the numerical procedure for 
solutions which had a saddle or a nodal critical point. 

We found that for a = 0.01 and 0.01 < a < 0.14 there are unique solutions for the 
sub-sonic, steady state flow with the effective speed of sound assumed to be constant. The 
critical point location varied monotonically in the range 2.8286 < x c < 3.1058 and the 
angular momentum constant varied monotonically in the range 1.00068 < b < 0.99528, 
i.e. we had x c ~ x ms = 3.0 and b 1.0. These results are not in any way affected by the 
supersonic flow at x < x c , i.e. no matter how complicated is the flow for x < x c it has no 
influence on the flow for x > x c . However, there might be a problem of the type presented 
by Abramowicz & Kato (1989) in their Fig. 2: the flow passing through the critical point 
may formally reverse at some radius x < x c , i.e. it may be globally impossible. This is not 
the case if the effective speed of sound is assumed to be constant also in the supersonic flow: 
our solutions continue smoothly all the way to the black hole, and the angular momentum 
remains almost constant within the supersonic flow, as expected. 

We were not able to find steady-state solutions for a > 0.14. Perhaps physical solutions 
cannot pass a nodal point in the slow direction, and the accretion is non steady for a > 0.14. 
We shall explore the consequences of variable effective speed of sound in section 5.2. 
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5. Some Analytic Results 

In order to have an analytic understanding of the behavior of the solutions and the 
dependence on the parameters a and a we attempt to replace eq. (14) by a simplified 
version of it, which retains the main topological features of the phase portrait, i.e. critical 
points and the 'Keplerian' curve, and also the dependence on a and a for small a's. First 
let us define v = v/a and u Q = v /a, where v , defined in eq.(18), is the 'Keplerian' radial 
velocity. Next we make the following assumptions 

x ~ 3, ^o — l, uj ~ 1, (19) 

which all follow from the assumption of a <C 1 and restricting the study to the vicinity of 
the 'Keplerian' curve (that includes the critical points). After these substitutions, being 
careful not to drop terms that are crucial for the main topological features, we end up with 

du ^ a(u-u ) 
dx VQau (u-l) 

We see that the 'Keplerian' curve, for which dv/dx vanishes, and the critical points, where 
v — vq — 1, are preserved. 

Next, we try to find a workable approximation for is (x). To do so, we notice that 
the interesting transitions in the nature of the physical solution (Fig. 6) happens where 
the critical points are close to the maximum of vq(x), which is also close to 1. With these 
approximations, we can use eq. (18) to get 

d\nu du 

^T"^ 1/3 - 2B(a; - 3) ' (21) 

where 

B- 1 = 32VEaa, (22) 

which can be integrated to give 

uq(x) ~ 1 + B[(x c - x b f -{x- x b f], x b = 3 + B"76. (23) 

x c is the position of the critical point, which is related to b by setting uj = 1 in eq.(15), and 
x b is the Saddle-Node bifurcation point, the boundary of the Saddle region and the Nodal 
region in Fig. 6. 

Now, let us define the new variables £, A and ji 

S = BW(x - x b ), n = B^VQa/a = [^- 3 } 1/2 , A = u - u . (24) 

62a 6 
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In terms of these new variables, eq.(20) becomes 
dA [j,- 1 A 



+ 2f . (25) 



Now, it is easy to find the eigenmodes M of the critical point, by setting A equal to 
M(£ — £ c ) and requiring that it satisfies eq. (25) to the first order in £ — £ c . This yields 

1 



M = - — {I - 4/i£ c ±yjl- 8fj£ c ), (26) 

which already gives the three saddle, nodal and spiral regions in Fig. 6. 

We do not intend to elaborate any further on the behavior of the solutions of eq.(25). 
This is mainly because the main property of the physical solution which is connecting 
the asymptotically 'Keplerian' curve to the critical point, is a global property and so its 
analytical investigation is not straightforward. The only clear conclusion that we can draw 
from eq.(25) is that the value of £ c for which the physical solution exists is a function of \x. 
The transitions in the nature of the critical point, points B (SN) and C (NN) in Fig. 6, 
that were discussed in the last section, occur at specific values of \x which can be already 
fixed from the numerical results in Fig. 6 and the definition of \x in eq. (24) 

Hsn ^ 1-22, fi NN ~ 0.53. (27) 

Of course the advantage is that now we can apply this result to all small values of a and 
not only a = 0.01. This yields 

^0~oT^ < ''014'' 3 ' No dal Point Through Slow Direction, 

ex. d ex 

3 < (777T7) < (t77^) 3 > Nodal Point Through Fast Direction, 



0.14 7 v 0.01 7 v 0.08' 



01 " 3 < (777rr)> Saddle Point. (28) 



v 0.08 7 v 0.01 

The numerical values calculated for a = 0.02 agree very well with these scalings. 



6. Discussion 

After numerical and analytical study of the properties of the solutions of eq. (14) in 
the last two sections, we are now ready to discuss the physical picture within the framework 
of our approximations. Of course the question of accuracy of these approximations should 
be addressed as well, and we intend to do so, at least in part. 
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Our numerical solutions fully confirm our guesses that the angular momentum constant 
~ L (i- e - b ~ f), and that dlogQ/dlogr pz —2 near r ms and for r < r ms (cf. eq. 9, and 
the text between eq. 8 and 9). Note: we have not assumed that I ~ l ms , this was obtained 
numerically as a consequence of two conditions: the disk had to be nearly 'Keplerian' at 
large radii, and the transition through the sonic point had to be smooth (cf. eqs. 16). 

6.1. Torque at the Sonic Point 

The starting point of this investigation was a recent controversy about the significance 
of the torque at the sonic point of a thin accretion disk. The traditional picture of Novikov 
& Thorne (1973) argues for a 'zero torque' boundary condition at the marginally stable 
orbit r ms . Historically, the possibility of a significant torque at the sonic point due to strong 
magnetic fields was first brought up in Page & Thorne (1974). This issue was pursued more 
seriously by Krolik (1999), Gammie(1999), Agol & Krolik (2000). 

Our disk model, with a smooth flow passage through the effective sonic point, has a 
small but non-zero torque there. The most relevant way to quantify the importance of this 
effect is to look at the energy generation rate due to differential rotation, which follows 
from the conservation laws: 




(cf. eq. 7). In the classical model of a geometrically thin disk with a no torque at r rns and 
'Keplerian' rotation for r > r ms , eq. (29) may be integrated to obtain 

^ = ££*= HO 55 = (-*)<*>>• £ °4' < 30 > 

Eq. (29) may be evaluated numerically for our disk model, and the results are shown 
in Fig. 13 for a = 0.01. A thick line corresponds to the classical model and two thin lines 
correspond to our models with a = 0.01 and a — 0.1, respectively. The differences are 
small but noticeable. The total energy released is reduced by fraction of one percent in 
the disk with a = 0.01 because the angular momentum constant is slightly larger than l ms , 
with b = 1.00068. The energy is increased by 2.5% for a — 0.1 model, because the angular 
momentum constant is slightly reduced, with b = 0.99636. 

Next, we calculated the change of accretion efficiency: 

Ae = L d 
eo (-M) c 2 e 



(31) 
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where L d was calculated with eq. (29) for 12 disk models with all combinations of the two 
parameters: a = 0.01, 0.02, 0.04, and a = 0.01, 0.05, 0.10, 0.14. The results are presented 
in Fig. 14 as a function of angular momentum constant: b = la/l ms - The important 
result is that the efficiency of accretion increases with increasing disk thickness, which is 
proportional to the parameter a, and with increasing parameter a. Within the range of 
parameters presented in Fig. 14 the efficiency is increased by up to 11% of e , i.e. from 

0. 0625Mc 2 to 0.0694Mc 2 . 

There are two related effects contributing to a small change in accretion efficiency. The 
fact that angular momentum constant Iq is not exactly equal l ms affects energy generation 
at all radii, even in the 'Keplerian' outer disk, as it is apparent in eq. (29). In addition, the 
disk near r ms is not completely 'dark' as it would be if the torque at r ms were exactly zero. 

The disk scale height H at the r ms may be estimated as 

(«) „(J?L) = (l) V \ xl . 63 a, (32) 
V r J ms \v ro tJ ms \3/ 

1. e. a = 0.04 corresponds to (H/r) ms 0.065. As we assume the effective speed of sound to 
be constant, the disk thickness increases with radius, and we reach H/r = 1 at r « r g /(2a 2 ). 
The 'Keplerian' condition at 'large radii' makes sense only for o <C 1, and this limits the 
range of parameters for which our disk models make sense. 



6.2. Is Nodal Passage in the Slow Direction Physical? 

Inspection of Figs. 2, 3, 4, 7, 8, 9, and 10 shows that the geometry of transonic flow 
appears to be similar when the critical point is of a saddle type or a nodal type, provided 
the solution passes the latter in the fast direction. The transition from the saddle to the 
nodal (fast) geometry is smooth, with no apparent problem. The sub-sonic flow connects 
well to a 'Keplerian' disk at large radii, and the supersonic flow approaches a free fall at 
small radii. This trouble free region corresponds to the aBC segment (thick line) in Fig. 
6, which ends at point C, where it contacts the borderline of the fixed points of the spiral 
type, and where a = 0.14. Obviously, no physical solution can pass through a spiral point. 
Above point C the analytical solutions pass the nodal critical point in the slow direction, 
and these solutions do not appear to be physical (cf. Fig. 11), as there is a sharp change in 
the slope of the solution (see the end of Sec. 4). We do not know if they are unstable, or 
perhaps no truly steady-state flow is possible for a > 0.14 (for a = 0.01). 

The study of the stability of these points with the assumption of constant v s and a 
is even less satisfactory than the steady state solutions. Nevertheless, Kato, Honma, and 
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Matsumoto (1988) show that, within these assumptions, local instabilities occur when 

a ss tt(r c ) > |^| c , (33) 

where ass is the Shakura-Sunyaev value of a, (ass — 2a; a). The dotted portion of the 
the physical curve in Fig. 7, which covers almost all of the sector Cd, shows the critical 
points that satisfy this criterion. This indicates that the nodal passage in the slow direction 
is probably unstable. We should point out that, contrary to the conventional picture 
(originally suggested by Matsumoto et al. 1984), this is not a generic property of the nodal 
critical points, but a generic property of the passage in the slow direction of the nodal 
point, where \dv/dr\ c in eq. (33) is small. 

We cannot prove that the transition through a nodal-type point in the fast direction 
is the only physically acceptable solution, i.e. that it is unique, but this seems to be likely 
upon inspection of Figs. 9, 10, 11, 12, and it also agrees with the finding of Artemova et 
al. (2001). If correct, this could resolve the ambiguity discovered by Matsumoto et al. 
(1984). Unfortunately, we cannot offer a simple physical explanation for the transition from 
saddle-type to nodal-type critical points while the parameter a increases. We stress that 
the physical nature of nodal points is outside the scope of our paper, which is concentrated 
on the issue of the inner boundary condition for geometrically thin steady state accretion 
disks. 



6.3. Discussion of Gammie (1999) paper 

There remains the issue which gave rise to this paper: why did Krolik (1999), Gammie 
(1999) and Agol & Krolik (2000) claim that the torque at inner boundary may be large 
even for thin, steady state disks, which remain thin even in the 'plunging region' ? Let us 
consider the case presented by Gammie (1999), as he provides the fullest model calculation. 
Gammie assumed that the disk is thin, c 2 s /c 2 <C 1...', and that a < 1 so that 
magnetic fields make a negligible contribution to the hydrostatic equilibrium of the disk 
He also writes: 'This picture leads one to consider a steady, axisymmetric inflow close to the 
equatorial plane of the Kerr metric' These assumptions are practically identical to ours. 

Some of the assumptions we made in this paper, and some made by Gammie are 
similar: we all assume that geometrically thin steady state accretion disks exist, and that 
magnetic fields are confined to the disk. Our model, with a constant and small effective 
speed of sound, implies that the flow remains geometrically thin also in the supersonic 
region, for r < r c . Gammie also assumed that the flow was geometrically thin everywhere. 
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In fact his flow thickness, H, was decreasing for r < r ms according to H/r = const <C 1 (cf. 
his eq. 4). 

A disagreement appears in the treatment of the transition from the 'disk' to the 
'plunging region'. Gammie adopts one set of assumptions for the flow in the 'disk', for 
r > r in , and a very different set of assumptions for the flow in the 'plunging region', for 
r < r in . At r in the model undergoes a dramatic jump in the adopted physical conditions: it 
is gas dominated and has negligible magnetic field for r > r in , and it has no gas pressure 
and it is magnetic field dominated for r < rj„. No justification is offered for this jump in 
physical conditions. The location of r in is adopted to be close to the marginally stable 
orbit, Tj n ~ "r ms . 

We do not introduce any rj„, and we treat the 'disk' and the 'plunging region' with 
the same equations. We require a smooth transition from sub-sonic flow in the 'disk' to 
supersonic flow in the 'plunging region' through a critical (sonic) point at some r c . We make 
no assumption about the value of r c , the condition of a smooth passage determines that 
value uniquely. We find that the critical point is located close to the marginally stable orbit, 
r c pa r ms , but it is in general not at r ms . The requirement of a smooth passage through the 
sonic point imposes a stringent restriction on the solution, it acts like a boundary condition. 

Another problem with Gammie's model is the assumption that the flow thickness is 
proportional to radius, i.e. that the flow becomes geometrically thinner in the 'plunging 
region', while the effective speed of sound rapidly increases with the decreasing radius in his 
model. This is incompatible with the hydrostatic equilibrium in the direction perpendicular 
to the equatorial plane. In fact there is no reference to the hydrostatic equilibrium in his 
paper. That equilibrium should hold as long as radial velocity remains much smaller than 
rotational velocity, i.e. as long as the accretion time scale remains much longer than the 
dynamical time scale. As the effective speed of sound increases in Gammie's model, the flow 
thickness should also increase, not decrease. A quasi-static expansion needed to maintain 
hydrostatic equilibrium is likely to conserve magnetic flux, i.e. B r H and B V H, where H 
is the disk thickness. As H increases the magnetic field gets weaker, B ~ l/H, and the 
magnetic torque at a given radius r is reduced: ~ Bj-B^rH ~ r/H. 

For the same reason the Alfven speed v a is reduced by flow expansion, so as to maintain 
approximate hydrostatic equilibrium: VA/v rot ~ H/r. 

It is beyond the scope of this discussion to speculate how the conclusions of AKG 
about the torque at rj n would change if the two major inconsistencies in Gammie's analysis 
were corrected: the jump at r in was replaced with full continuity of all physical quantities, 
and the consequences of hydrostatic equilibrium in the 'vertical' direction were incorporated 
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in the model. These two problems are serious enough to doubt the conclusions: the torque 
at r ms may be large and the accretion efficiency may exceed 100%. 

No direct comparison is possible between our model and Gammie's (1999) model, as the 
latter does not specify the value of H/r. Amazingly, his results appear to be independent of 
H/r value, while our results are sensitive to the assumed disk thickness, as given with the a 
parameter (cf. eq. 32, and Fig. 14). 

The consequences of ad hoc assumptions may be very diverse. AGK claim the 
importance of magnetic fields increases with decreasing radius, and becomes dominant in 
the 'plunging region'. Using different ad hoc assumptions Li (2002) comes to the opposite 
conclusion in his simple analytic model. 

A large part of Agol & Krolik (2000) analysis remains correct: there may be a torque 
applied to the inner edge of the accretion disk provided there is a large scale magnetic 
field, possibly threading the black hole. This is a modified Blandford-Znajek mechanism 
(e.g. Li 2000, Wang et al. 2002, and references therein). The gaseous disk may remain 
geometrically thin while the stresses are transmitted by a large scale magnetic field which 
has a vertical scale height comparable to radius. However, this picture is conceptually very 
different from the case considered by us and by Gammie (1999), as we and Gammie assume 
that the magnetic fields are confined to a geometrically thin flow. 

7. Conclusions 

The condition given with the eq. (1) of this paper is local, and it is approximately valid 
even if the parameter a varies with radius. It requires the assumption of a steady state to 
hold. However, we think it is reasonable to expect that it may also hold for a disk with an 
accretion which is steady in a statistical sense only, as long as the disk thickness does not 
exceed the value H <C r throughout the fluctuations cycle. Obviously, it is important to 
verify this claim with the full 3-D time dependent simulations. However, current numerical 
models are limited by the available computer power to disks which are much thicker than 
the one considered in this paper. They are also limited by the absence of cooling processes 
which are essential for the formation of thin disks. At this time there is no observational or 
theoretical reason to exclude the possibility that geometrically thin disks exist. Currently, 
there is also no way to prove their existence, as the relevant thickness is that of the layer 
over which the magnetic fields transferring angular momentum extend, and this quantity is 
not readily observable. 

We have found a unique steady state solution for a simple model with a constant 
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effective speed of sound v s = 0.01c, and a given viscosity parameter a < 0.14. The solution 
was found by approaching the critical point from large radii, i.e. from the sub-sonic flow 
side, where at large radii the disk was assumed to be 'Keplerian'. The same solution 
was found starting integrations from the critical point and proceeding out and seeking a 
solution approaching 'Keplerian' at large radii. No assumption was made about the value 
of the angular momentum constant Iq, which was determined by our choice of boundary 
conditions: 'Keplerian' disk at large radii, and a smooth passage through the sonic point. 

Our numerical model confirms the validity of our eq. (1): we find that in a thin steady 
state disk the sonic point r c is located close to the marginally stable orbit r ms , and the 
value of angular momentum constant Iq is close to l ms . 

The supersonic flow has no effect on the critical point, as expected on general grounds. 
In particular, we demonstrated that, for geometrically thin inflow, with constant effective 
speed of sound, there is no global problem of the type envisioned by Abramowicz & Kato 
(1989) in their Fig. 2. Once the flow passes the critical point it plunges into the black hole 
along a spiral, roughly conserving angular momentum, as asserted in many papers written 
in the 1980s. 

In our disk model the torque at the sonic point is small, but not exactly zero, and 
the total disk luminosity is modified, but only by several per cent. This small torque 
and a modest change of accretion efficiency makes the largest difference near r ms , where 
the accretion flow is no longer 'dark', as it was under strict 'no torque' inner boundary 
condition. This effect should be included in calculating spectra of thin disks. 

We emphasize that for a disk to be thin its thickness H must include the magnetic 
structures responsible for angular momentum transfer, not just the gas layer. The effective 
speed of sound v s is affected by the total pressure, which includes the contribution of 
magnetic pressure (cf. eq. 5b). 

Just as discovered by Matsumoto et al. (1984) we also found that the critical point 
is of the saddle-type for small values of a, and becomes a nodal-type for large a, in our 
case for a > 0.08. Unfortunately, we cannot provide a simple physical explanation for 
this transition. However, we presented plausible arguments that the transition through 
the nodal-type critical point is unique, as there is only one solution which appears to be 
smooth. This was also demonstrated by Artemova et al. (2001). 

We found that our steady state solutions exist only for a < 0.14(100f s /c) 1 / 3 , and it 
appears that there are no physically sensible solutions for a > a cr u ~ 0.14 (lOOfs/c) 1 / 3 . 
It is not clear if this is a general result or just an artifact due to simplicity of our model. 
The issue is not the particular value of a crit , which is certainly model dependent, but the 
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very existence of a crit . Unfortunately, we cannot offer a simple physical explanation for this 
finding. 

Krolik (1999), Gammie (1999) and Agol & Krolik (2000) claimed that their 
geometrically thin steady state disks had large torques at r ms , and that accretion efficiency 
could be larger than 100% of Mc 2 , i.e. a thin accretion flow could extract energy from a 
spinning black hole. This may be correct if there are large scale magnetic fields, like those 
proposed for the modified Blandford-Znajek mechanism, also referred to as a magnetic 
coupling model (e.g. Agol & Krolik 2000, Li 2000, Wang et al. 2002, and references 
therein). However, if the magnetic fields are confined to a thin accretion flow, as in the 
model described by Gammie (1999), the very high accretion efficiency is almost certainly a 
consequence of several inconsistent assumptions (cf. our Sec. 6.3). It would be interesting 
to learn how Gammie's results change when his errors are corrected, but this is beyond the 
scope of this paper. 

With our thin disk model we found that there is a small torque at r ms , and there is 
a corresponding modest change of accretion efficiency. It is not possible to compare our 
results directly with those obtained by Gammie as he does not specify how thick his disk is. 

We are very grateful to the anonymous referee who relentlessly pressured us to refine 
our analysis. It is a great pleasure to acknowledge many useful and critical comments by 
Dr. J. Goodman. 
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A. 'Quasi-adiabatic' accretion 

Our assumption that the accretion flow remains geometrically thin at all radii: in the 
sub-sonic part for r > r cr , and in the supersonic part for r < r cr , is physically equivalent to 
the assumption that magnetic fields dissipate effectively at all radii. In this Appendix we 
discuss the consequences of a possibility that there is no field dissipation at radii smaller 
than some transition radius r tr . In such a case the magnetic field pressure and the disk 
thickness may increase for r < r tr - Krolik (1999) pointed out the possibility that the 
differential rotation within the supersonic flow may increase the effective speed of sound so 
much that we expect the Alfven speed in the fluid frame to be ~ c ...' (Krolik 1999). 
Obviously, if the magnetic field dominates the flow then the effective speed of sound is the 
same as Alfven speed, and it also becomes relativistic. 

The largest energetically possible increase of the effective speed of sound may be 
obtained by assuming a 'quasi-adiabatic' flow: all energy transfered from differential 
rotation into the accretion flow is used to increase magnetic energy, no magnetic energy 
is dissipated, and none is radiated away. This is certainly a limit which is not likely to 
be reached in realistic flows. Nevertheless, we explored this possibility as a limiting case, 
and we refer to it as a 'quasi-adiabatic' approximation. Note that the model developed by 
Gammie (1999) was also 'quasi-adiabatic'. 

We repeated the calculations of our models with a = 0.01. In the first case we adopted 
a = 0.01, in the second a = 0.1, and we obtained: (x c = 2.82859, b = 1.00068) and 
(x c = 3.04505, b = 0.99636), for the two cases, respectively. Obviously, these values were 
identical to the models described in Section 4. We continued these solutions some distance 
into the supersonic part of the flow, down to a transition radius x tr , where we abruptly 
switched the flow from having constant effective speed of sound v s = 0.01c, and constant a, 
to a 'quasi-adiabatic' approximation. All physical quantities varied continuously, but their 
derivatives experienced a jump. 

We assumed that at x tr the ratio of magnetic pressure to total pressure was equal to 
a. For r < r tr the energy of differential rotation was pumped into tangled magnetic fields, 
hence the magnetic pressure increased. No magnetic energy was dissipated, and the gas 
pressure was assumed to change adiabatically, i.e. the gas pressure decreased because of 
decompression. Therefore, the ratio P m ag/Ptot increased down the flow for x < Xt r , and 
the parameter a increased correspondingly. We assumed that hydrostatic equilibrium was 
maintained in the direction perpendicular to the disk equatorial plane, as that equilibrium 
was established on a fraction of rotational period, while the flow followed many rotations 
prior to crossing the black hole horizon. The gas density rapidly decreased as the flow 
accelerated toward the black hole and expanded in vertical direction. The vertical expansion 
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of the flow was a consequence of the increase in the effective speed of sound, which became 
practically equal to the Alfven speed. 

We found that if the transition from v s = ac = const to a 'quasi-adiabatic' flow was 
made sufficiently far into the supersonic region, with x tr < 2.79 and x tr < 2.78 for the 
two models, respectively, the flow proceeded smoothly all the way to the black hole, and 
remained similar to the flow in our original model. At the transition radius, x tr ~ 2.8, the 
flow was moderately supersonic, with the effective Mach number of 1.3 and 2.4, for the two 
cases, respectively. Down the flow, for x < x tr , the magnetic pressure increased, and the 
parameter a also increased to a ~ 1. The effective speed of sound gradually increased by a 
factor ~ 3, increasing the flow thickness correspondingly, but the flow remained supersonic. 

When the transition was made at a larger value of Xt r , i.e. at a smaller effective Mach 
number, the flow pattern was initially 'normal', i.e. the a, the effective Mach number 
and the Alfven speed all increased gradually with decreasing radius. However, at still 
smaller radii the flow first became sub-sonic, and later the numerical solution 'reversed' 
the flow direction, changing the sign of v r , i.e. the flow proceeded back toward large radii. 
Obviously, this is physically impossible, as pointed out by Abramowicz & Kato (1989). 
Superficially a global steady state solution could not exist for the transition radius larger 
than approximately 2.8r 9 . However, it is virtually certain that before the flow would 
'reverse', the deceleration would generate a shock, and a transition to a sub-sonic accretion. 
Perhaps a steady state flow could be maintained, with the second sonic point located down 
the flow from the shock location, and the matter ultimately plunging into the black hole. 
However, these complexities are beyond the scope of our paper. 

The results of our numerical experiment are easy to understand qualitatively. As 
the dissipation of magnetic energy is switched off ('quasi-adiabatic' approximation) 
the magnetic pressure builds up, and the flow expands in vertical direction. Had the 
transition been made at a large radius, the disk would become geometrically thick. The 
very assumption that the disk is thin is obviously equivalent to the assumption that the 
dissipation of magnetic field is efficient. 

Gammie (1999) assumed that there was no magnetic dissipation in his model i.e. his 
model was also 'quasi-adiabatic'. But there were major differences between his and our 
models. Gammie assumed that H/r = const <C 1 at all r < r in ps r ms , while our flows 
thickened appreciably for r < r tr < r ms . Gammie assumed that magnetic fields had a 
very special smooth spiral geometry, while we assumed they were tangled. His solutions 
for the flow were very different from ours. The main result of our experiment is that 
in a 'quasi-adiabatic' approximation the flow thickens, which is hardly surprising. We 
also found that properties of a 'quasi-adiabatic' flow strongly depend on the value of the 
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transition radius r tr = x tr r g . Gammie did not explore this diversity, effectively adopting 

fin ^tr Tms- 

We were curious to study the consequences of a 'quasi-adiabatic' approximation, and 
we were intrigued to find that even within this approximation our solutions were very 
different from those obtained by Gammie. However, this experiment is not relevant to the 
main topic of our paper, which was focused on the properties of thin accretion disk. 
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Fig. I. — The angular momentum constant b is shown as a function of the position of the 
critical point x c for a — > 0. We see that, for b > 1, there are two fixed points, at b = 1, there 
is one degenerate fixed point at x = 3 (marginally stable orbit), and for b < 1 no fixed point 
exist. 




Fig. 3. — The Phase Portrait for a 
at x — 2.94 and a nodal point at x 



= 0.01, a = 0.15 and b = 
= 3.1. 



0.9949. There is a saddle point 




Fig. 5. — The Phase Portrait for a 
for these parameters. 



= 0.01, a = 0.15 and b 



= 0.9927. There is no fixed point 
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Fig. 6. — The parameter space of the problem for a = v s /c = 0.01. x c is the radius of 
the critical point in units of r g and a is the viscosity parameter. The curve aBCd shows 
the position of the fixed points that allow physical solutions (cf. Sec. 4), while the dotted 
section of the curve are unstable solutions according to eq.(30). Also the parameters of phase 
portraits in the other figures are represented by numbered diamonds. The number of each 
diamond is the number of the associated figure in this paper. The Phase Portrait presented 
in Fig. 5 has no fixed points, and therefore it is not represented in this figure. 



Fig. 9. — The Phase Portrait for a = 0.01, a = 0.1 and x c ~ 3.047. The physical solution 
(thick curve) passes a critical point of the nodal type, in the fast direction. 




X 



Fig. 10. — The vicinity of the critical points in the phase portrait of Fig. (9). 




Fig. 12. — The vicinity of the critical points in the phase portrait of Fig. (11). 
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Fig. 13. — The energy generation rate in the classical model is shown with a thick solid line, 
and for our a = 0.01 disk model with two thin lines, corresponding to a — 0.01 and a — 0.1, 
respectively. The disk luminosity is normalized to unity for the classical model. 
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Fig. 14. — The relative change in the accretion efficiency for our disk models with constant 
effective speed of sound: v s = ac are shown for three values of a (0.01, 0.02, 0.04) and for 
four values of the parameter a (0.01, 0.05, 0.10, 0.14). For a = 0.01 the angular momentum 
constant Zo is smaller than l rns and the accretion efficiency is reduced compared to the model 
with exactly zero torque at r ms . The efficiency of accretion increases with increasing a and 
with increasing disk thickness, which is proportional to a. 



